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I. INTRODUCTION 

The pair plasmas consisting of only positive- and negative- charged particles of equal 
mass have attracted special attention mainly because of the astrophysical applications. In 
the early universe during the lepton era, ultra-relativistic electron-positron (e-p) pairs con- 
tribute largely to the matter contents of the Universe jl| . The gamma-ray bursts - the most 
concentrated electromagnetic explosions in the Universe - are believed to be related with 
the enormous energy release in compact regions on short time-scales. This energy release 
leads to the formation of a highly dense optically thick e-p plasma that expands and cools 
down remaining relativistic I2I. Such pair plasmas exist also in active galactic nuclei, in the 
relativistic jets, j^, and in the pulsar magnetospheres 4]. ^ 
Although there are many laboratory setups in which e-p pair plasmas are produced [5!|- 







one of the most interesting recent laboratory accomplishments is the successful creation 
of " sufficiently" dense pair-ion (pi) plasmas - the first such plasma consisted of equal-mass, 
positive and negative fuUerene ions {Cqq and Cqq) j^. 

Unlike the e-p plasma systems (both of the astrophysical and laboratory variety), the 
fuUerene plasma has a long enough life time that the collective behavior peculiar to the 
plasma state can be experimentally investigated under controlled conditions. The frequencies 
associated with the collective modes (plasma frequency, acoustic and Alfven frequencies) in 
such plasmas tend to be rather low. Fortunately, the group of Hatekayama and Oohara have 
already made considerable progress in the production of the hydrogen, H~^- H plasmas 



id . . Since the initial report, both the quality and quantity of this light pair ion plasma 
has been steadily improving [12]. 

Since many properties of pair plasmas (a symmetric pair plasma, for instance, cannot sus- 
tain charge separation) are different from the ordinary electron-ion (e-i) plasma, controlled 
experiments would not only advance fundamental physics, but also create a laboratory to 
simulate and understand a variety of phenomena taking place in astrophysical environments. 
A basic requirement for long-time-scale experiments, will be that the pair annihilation time 
scale is many orders of magnitude larger than the plasma period. 

Theoretical investigation of pair plasmas have followed two distinct tracks: 
1) the first track emphasizes the special properties that stem from the symmetric (pair 
particles have the same inertia, temperature etc.) nature of the pair plasma. This is a highly 



studied field Il3l 



both in the astrophysical as well as the laboratory contexts. In 



the e-p plasmas, problems like solitary structure formation are studied for instance in 



while for the ion pair plasmas much of the linear as well as nonlinear work has been devoted 



to understanding and interpreting experimental results (see e.g. |18l . Il9|). in particular the 
dispersion curves, described in [91] 

2) The second group of papers deals with a whole lot of interesting phenomena that arise 
because the symmetry of the pair plasma is mildly broken through some mechanism which 
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creates some disparity between the constituent fluids 
of the more interesting results of symmetry breaking is the creation of localized nonlinear 
structures; the particular properties of the structure will, naturally, depend upon the mode 
of symmetry breaking. 

Although it is natural to imagine that pair plasmas have to be symmetric: that the 
charges, densities, temperatures and masses of positively and negatively charged particles 
are equal. But both observations and experiments indicate that asymmetry may appear, 
quite naturally, at some stage of their evolution. In fact asymmetry could also be engineered 
in experiments if such a state could show interesting properties. Different species, not pro- 
duced in identical conditions [l^, for instance, could have different thermal speeds (tem- 
peratures). One could also arrange experiments with different setups for different species 
when, for instance, there are fractions of heavier/lighter ions or there is a mixture of different 
mass/temperature species with opposite charges. This way one could mimic the conditions 
pertinent to astrophysical pair plasmas. 

A much investigated example of broken symmetry is a pair plasma contaminated by a 
small fraction of charged particles with different mass (lighter/heavier than main species). 
Symmetry breaking could also occur when the constituent elements of the two fluids have 
slightly different masses, or the fluids have slightly different temperatures. As expected, 
symmetry breaking induces new properties that are different from what pertains for pure 



pair plasmas. It was shown in 
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26l | that pair plasma, contaminated by the heavier 



immobile ions, can support 3D stable completely localized structures of EM radiation - 
"light- bullets", non-diffracting and nondisspersive EM pulses of pancake shape with large 
density bunching. Existence of such localized structures is not possible in pure pair plasmas. 
Localized nonlinear structures of EM radiation were also found in a hot e-p relativistic 
plasma containing a small fraction of cold electron — ion component [27!]. Similar behaviour 
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could be expected in doped (or dust-contaminated) fullerene plasmas in laboratory [28 1. 

The present study concentrates on establishing the existence of electromagnetic (EM) 
solitonic structures in pair plasmas that could have been created due to asymmetries of 
different origin. We will work out the consequences of two classes of symmetry breaking; 
asymmetry arising from a small temperature difference of the constituent species and the 
asymmetry arising from a small difference in the species masses. The plasma is assumed to 
be underdense while the EM pulse is longer then characteristic skin length of the plasma. 
In electron-ion underdense plasma the formation of solitonic structures takes place only 
at the frequency close to the plasma frequency since in such plasma Raman instabilities 
dominate the process of soliton formation. In contrast, in pure symmetric pair plasma, 
ponderomotive forces are same for different species and the excitation of longitudinal waves 
by the EM pulse and Raman instabilities can not develop. In slightly asymmetric pair 
plasmas the ponderomotive forces acting on positively and negatively charged species are 
slightly different; the generation of weak ambipolar electrostatic potential takes place and the 
effects related to Raman instabilities can be ignored to leading order. This potential plays 
a fundamental role in structure formation; it acts as the binding "glue" that concentrates 
matter and radiation in a small region. The nonlinearity due to temperature asymmetry is 
found to have a focussing-defocussing form (different from the nonlinearities originating in 
other modes of symmetry breaking) imparting a rich structure to the corresponding solitons. 

II. MODEL 



Let us assume that the velocity distribution of particles is locally a relativistic Maxwellian. 
Then the dynamics of the fluid of species a [a represents negative and positive particles of 
any origin) is contained in the equations (see for details |29[|): 

Q 

Q^{GaPa) + "^oaC^V(G'„7«) = e^^E + (u« X n^), (1) 

whose curl converts the equation of motion to the vortex-dynamical form: 

dfl 

= V X (u« X n„), (2) 

where = 'ja'tTT-a'^a is the hydrodynamic momentum, E and B are the electric and 
magnetic fields and fla = (ea/c)B + V x GaPa is the so called generalized vorticity. 
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Here Uq, denotes he hydrodynamic velocity, 7q, = (1 — u^/c^)~^/^ = (1 + p^/mg^c^)^/^ 
is the relativistic factor and rriaGaiza) = iTia K^{za) / K2{za) is the thermally enhanced 
"effective mass" , [2;^ = ""^aC^/^a] 5 where K^, are the modified Bessel functions and rria and 
Tq, are the particle rest mass and temperature of species a, respectively. For nonrelativistic 
temperatures (T^ <^ m^c^) Ga = l+5Ta/2maC^ and for the ultra-relativistic temperatures 
(Tq 3> rriaC^) Ga = ^Ta/ma(? ^ 1. Note, that the relativistic thermal pressure Pa 
[= {na/'^a)Ta in the rest frame with being the density in the laboratory frame of the 
pair-fluid] appears through the temperature dependent factor Ga defined by 7q:VPq = 
ma(?noSGa- The system of Eqs. (II])-® is augmented by the equation of state: 

— , exp[-z„K2(za)] = consta, (3) 

which yields the usual results {n'^T^^'^ = const in non-relativistic limit for mono-atomic gas 
and n'^T^ = const in ultrarelativistic case for photons). 

From ([2]) it follows, that if the generalized vorticity is initially zero [fla = 0) everywhere in 
space, it will remain zero for all subsequent times. We assume that before the EM radiation 
is "switched on" the generalized vorticity of the system is zero. 

And for both species we have the continuity equation: 

^ + V ■ (naUa) = 0. (4) 

To study the nonlinear propagation of intense EM wave in a pair plasma (could be 
relativistically hot) consisting of negative and positive ions we must couple the equations of 
motion with Maxwell equations. In terms of the vector (A) and electrostatic (0) potentials 
defined by: 

1 (9A 

E = — — -V^; B = VxA, (5) 
c ot 

the basic equations take the form (Coulomb gauge V ■ A = 0): 

d^A ... d 



A A + c- ( V^) - AircJ = 0, (6) 
and 

Alp = -47rp, (7) 
where for the charge and current densities we have respectively: 

P = ^<^ana; J = XleanoUo. (8) 
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Our purpose is to explore the possibility of finding localized structures in such plasmas; 
the equilibrium state is characterized by charge neutrality: Uq — Uq, where Uq and Uq 
are the positive and negative charge unperturbed densities. The subscript a hereafter will 
indicate the negative (a = — ) and the positive {a = +) ions. In terms of dimensionless 
variables: 
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t — UJ-t , 



(9) 



where cv- — {Aire'^nQ /m )^/^ is the Langmuir frequency of negative species, the entire set of 
defining equations reads: 



^-AA+-(V0) + 



r+ 



0, 



+ vr± 

at 



dA 



dt 



0, 



(10) 

(11) 

(12) 
(13) 



where it was convinient to introduce temperature dependent momentum 11='= = C^p"^ and 
relativistic factor r± = G^'y^ = ^ {G^f + (n±)2. The equations of state for species then 
read: 

= ^ (14) 

with /(r±) = [r±X2(l/7'^)/G'^]exp[G±/T±] and the equilibrium temperature of 
species. 

We now discuss the propagation (along the z axis) of a circularly polarized EM wave with 
a mean frequency Uo and a mean wave number ko 

1 



-(x + iy)A(r_L, t)exp{ikoZ — iujot) + c.c, 



(15) 



where A{z, t) is a slowly varying function of z and t and x and y are the standard unit vec- 
tors. The choice of circular polarization is dictated by considerations of simplicity (prevents 
harmonic generation). We shall now follow standard methods to analyze the system. We as- 
sume that the longitudinal extent of the pulse is much shorter than its transverse dimensions 



7 



{dA/dz > V_L^)- The gauge condition gives us = (V^o)(V_l ■ A_l); I^^I < \A±\. Con- 
sequently the effects related to A^ will turn out to be negligibly small. Then, in the slowly 
varying amplitude approximation, the transverse component of (fT2il is integrated yielding: 



n 



± 



(16) 



where the constant of integration is set equal to zero since particle hydrodynamic moments 
are assumed to be zero at infinity where the field vanishes. Note that for the longitudinal 
motion the equations of motion can be treated one-dimensionally. 

Longitudinal dynamics is described by the z component of the equation of motion (fT2|) 
and continuity equation. We note, that due to the circular polarization of EM wave 7^ = 
[l + \A\'^/{G^f + {Ilff/{G^fY/^ doesn't depend on the fast time uJq^ and all the quantities 
in ( |T2l) vary on a slow time scale. Therefore, we can introduce the following variables for 
convenience: ^ = z — Vgt ,t = t, where Vg = ko/ujQ is the group velocity of the EM wave 
packet. Assuming Vgd/d^ ^ d/dr, straightforward algebra gives the following integral of 
motion: 

r±\2i 1/2 



1 + 



Ul 



Vgllz ± 



const 



± 



{G^y ' (G±)2_ ^^."^^ • (17) 

The constants of integration are determined from the boundary conditions: the EM fields 
and plasma momenta vanish at infinity; they are found to be G^{T^). 

In this paper we will deal with transparent plasmas, i.e, c^o ^ 1, and Vg ^ 1; such plasmas 
are found both in astrophysical and laboratory conditions (non-transparent plasma case will 
be discussed in a future publicarion). From the continuity equation we have: 
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(18) 



and straightforward algebra leads to: 
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where G^ = G^{T^). Equation ( fTTl) yields: 
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which allow us to write densities fully in terms of potentials and A: 



n 



2Gt 



1 + 



It 



Gt 



1 



(22) 



(23) 



In pure pair plasma with equal temperature species, the radiation pressure gives equal 
longitudinal momenta to both the negative and positive ions (since their effective masses 
are equal [G" = G^ = G]) and thus does not produce any charge separation («_ = 
and = 0) (Berezhiani and Mahajan, 1994). Berezhiani and Mahajan (1995) showed that 
the situation changes by introducing a small fraction of heavy ions; with the "symmetry 
breaking" between hot electrons and positrons, it becomes possible to generate a finite 0. 
Several studies showed the existence of electrostatic potential due to the small fraction of 
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22| : pair plasmas respond similarly 



27|. 



different temperature electrons 

In 26| it was demonstrated that in the pair ion plasmas, for which the symmetry is 



broken by a slight contamination (doping) through a heavier immobile ion, the electrostatic 
potential is no longer zero and such plasmas can support stable localized EM wave struc- 
tures even in the nonrelativistic limit appropriate to the current and near future laboratory 
experiments. As we will see below, the creation of electrostatic potential is also possible due 
to the difference in initial temperatures for different (major constituent) species (yielding 
the difference in "effective masses" of species). 



III. FORMATION OF LOCALIZED STRUCTURES IN PAIR PLASMAS WITH 

TEMPERATURE ASYMMETRY 

We introduce the "asymmetry" through temperature difference between the two species. 
The existence of temperature asymmetry is both experimentally and observationally justified 
(see the references given in the Introduction) and leads to different "effective masses" even 
though the real masses are equal (m~ = m"*"). We would expect that this asymmetry in 
"effective masses" would lead to effects similar to the ones studied for ordinary e-i plasmas. 
As a first step we assume, that temperatures are only slightly different 

^ = ^ 1, (24) 
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f3 is close to 1 so that we can have small but finite 
This assumption lets Eqs. fl20|) to be written as 
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(26) 



where the following normalization g"^ = (G^/G^) and = (pjG^ were used. Eq. ([23 
written fully in terms of potentials 0, A, relativistic factors and /5, transforms to: 



n 



1 + 



n 



(1 - /30)2 
, {9-f + \A\' 



(27) 

:i+0)2 

where A = A/G'^. From now on we will omit the and it will be assumed that (j) and A 
are normalized to G^. 

For clarity of exposition we will not attempt to analyze the general case confining ourselves 
to the two extreme limits: non-relativistic and super-relativistic temperatures. 

To close Maxwell equations, we need to evaluate = n^n^/F^, the current that goes 
into ( ITOl) . Using Eqs. equations ( l25l) . ( I26l) and the normalized version of Eq. ( |T6l) 



the expression for the total current is easily derived to be: 



+ 



A 



(28) 



(29) 



.1-/30 ■ 1 + 

Using the quasi- neutrality condition = n~ (characteristic length-scale of wave L ^ 1) and 
plugging it into the equations (127|) . and taking into account that for a transparent plasma, 
heating of both fiuids is very weak (implying ~ 1) we can derive ~ (1 — /5)'0(|Ap) 
(where ^/'(lAp) < 1). Putting all the pieces together (evaluating the current, actually the 
(J^ — 2A) in our notation) the nonlinear term in Eq. (fTOj) is expressible as: 



NL 



+ 



1 



-0[(l-/?)-2/30], 



(30) 



1-/50 ■ 1 

where we do not neglect w.r.t. (1 — /5) ^ 1. To complete Maxwell's equations, we have 
to relate the electrostatic potential (< 1) with its source, in this case the temperature 
difference between the species. 
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A. Super-relativistic temperature pair plasmas 

For super-relativistic temperatures T"^ ^ 1 and = 4T='=(:» 1), g"^ = T'^/T^, and Eq. 
Ml) reads: 



9^'. (31) 



/ \ oo , 

which, in conjunction with ( l25l) . ( l26l) and ( 13T1) . yields : 

3^ = (1-/50)1/2 ' = (1 + 0)1/2 (32) 

evaluating the "effective mass" self-consistent ly in terms of 0. Straightforward algebra 
and the use of quasi-neutrality condition yields 



"[(l-/?)-2/3, 



leading, finally, to 

2 





\A\ 
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(l + /3/« 


\A\ 





with K = - . (33) 

We see that ~ (1 — /9) when |y4p ^ 1 and < (1 — /3) when |y4p < 1 and our estimation 
of nonlinear term (!30!) remains valid. 

Note, that with electrostatic potential defined by fl33l) one learns that that heating/cooling 
of both fluids is weak {g+ ~ (1 + 0.5/? 0); ~ (1 - 0.5 0)). Also from the Eqs.([27D the 
species densities come out to be: 

n±~l + ^|A|2. (34) 



B. Non-relativistic temperature pair plasmas 

For non-relativistic temperatures T^,T^ <^ 1 and = 1 + 5T^/2 . In this limit the 
relevant relations are: 

s-^l + l(T--T-). (35) 

and 

P = ^^l + l{T--T^). (36) 
Straightforward but tedious algebra leads us to 



g± = l + 



(37) 
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where 

ff+ = ^/3T+0 and /7_ = -^T-0. (38) 



From (EZD-dHHl) we find: 



2 ' oo r 2 



g+-g-^\(t>T:^{l + l3), (39) 



from wliicli 



^ 0T-(l + /3), (40) 

readily follows. Repeating the same procedure as we adopted earlier, we may derive 

(1-/3) kUP 1 

^ '^^ ' ' with K = - (41) 
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2 (l+/5fi:|A|2)' 2 

an expression that has the same general form as Eq. (133|) that pertains to the super-relativistic 
temperature regime! In fact, the two limits are very similar - for either case the plasma 
heating/cooling is weak and the density bunching ra^ ~ 1 + /5 |ylp/2. 



IV. LOCALIZED STRUCTURES 

Our interest, here, is to delineate the conditions for the existence of localized structures 
in pair plasmas. For the temperature asymmetry driven system, the perpendicular current 
is assembled from Eqs. ( !33l) and (HTj) : 



j^^-(2-(j)[{l-(3)-2(3<j)\) A 
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(42) 



with K = 1/2 for non-relativistic temperatures and k = 2/3 for relativistic temperatures. 
The final equation for the vector potential (equivalent to the final defining equation (43) of 



251]), derived from Eq. (fTOj) and the preceding expression of current, reads: 



dA 1 + /5 d^A _2 . 
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A = 0, 



(43) 



with 



;i-/9)'«i, 
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where we have redefined m_ as m_ — > m_G^ that introduces the "effective mass" for 
negatively charged ions. The wave frequency ujq satisfies the dispersion relation: uJq = 
+ {1 + (3) (in dimensional units this reads as: = + {1 + (3) uj'i ) . In Eq. (H3ll 
the weak dependence on the transverse coordinates has been retained. Note, that in spite 



of the fact that dA/d^ ^ V±A the second and t 
of the "transparent plasma" {uJq ^ 2) condition 



le third terms can be comparable because 



30|. 



With self-evident renormalization, Eq. fl43l) can be written as: 

^^ + ^ + ^lA + H\A\')-A = 0, (44) 
with the nonlinearity function given by: 

Equation (jH]) is nothing but the nonlinear Schrodinger equation (NSE) with a saturating 
nonlinearity. The saturation function F(|y4p) constitutes a "new" type; it has an unusual 
form in the sense that in the ultrarelativistic case |Ap ^ 1 it tends to vanish. To the best 
of our knowledge NSE with such a "vanishing saturation potential" has not been derived 
and reported so far for any known physical system. 

From the refractive index 611^1 = F{I) (/ = |y4p is the intensity of the EM field), we may 
deduce that the plasma is self-focusing {d{6nni)/dl > 0) provided / < 1 while it becomes 
defocusing {d{Snni) / dl < 0) for higher intensities (/ > 1). For a localized intense EM pulse 
with a peak intensity > 1, the medium, thus, reacts differently to different parts of the 
EM pulse - focusing the peak region while defocusing the wings. 

An immediate consequence of the fucusing-defocusing saturating nonlinearity (originating 
from the temperature-asymmetry driven mechanism) can be illustrated by considering a 
modulation instability of quasi-monochromatic EM wave. Equation (jH]) is satisfied by the 
plane wave solution A = exp(zrF(|y4oP)) + c.c. The standard stability analysis then 
shows that a linear modulation with frequency Q and wave number K obeys the dispersion 
relation fi^ = K^[K^ — 2^0(1 — Al)/{1 + A^)^] which exhibits a purely growing mode 



if < 1 and K < K^r = y2AQ{l — Aq)/{1 + AqY, while for the ultrarelativistic case 
{Aq > 1), there is no modulation instability. One can expect that the modulation instability 
of moderately intense field {Aq < 1) in the nonlinear stage will lead to the break up of the 
field into soliton-like pulses with a characteristic length corresponding to the optimum scale 
of instability (~ \/2/Kcr)- 
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In what follows we investigate the possibility of finding stable solitonic solutions of Eq. (H^ 
under a variety of conditions. For stationary solitons, we look for solutions that are "spher- 
ical" symmetric: A = A{r) exp(iA) where A is a constant measuring the nonlinear frequency 
shift. Notice that the co- moving coordinate ) can be treated on an equal footing with 
the spatial coordinate (r^). In terms of the radial variable r = (r^^ + ^^)^'^^, Eq. (H^ re- 
duces to an ordinary differential equation that can not be analytically solved. However, to 
better understand the results of possible simulations it is helpful to rewrite it (after trivial 
manipulations) as the equation describing a "particle" moving with friction in the potential: 

where the "effective potential" isV{A) = -XA^ + \n{l + A'^)-A'^/{l + A'^). Here D(= 1,2,3) 
dimension of the problem. 

The profile of the potential for different values of the nonlinear frequency shift A is 
presented in Fig.l: 1) curve "a" in this plot corresponds to A > A^J;^) ^ 0.2162, 2) A = A^J;^) 
for the curve "b", and 3) the potential in the range < A < A^J.^^ is given in curve "c". 
The solitary solutions correspond to the effective particle which at the "moment" r = 
rests at a point with coordinate Am- Then it rolls down (r > 0), dissipates "energy" and 
approaches asymptotically (r oo) the potential maximum at A = 0. It is obvious that 
solitary solution can not exist for A > X^r^^ while for A = A^J.^^ solitary solution does exist. 

However in the range covered in "c" (0 < A < A^J.^-') solitary wave solution exists in 
more than one dimension; the corresponding Am for a given eigenvalue A has to be found 
numerically. Notice, that in 1-dimension (1-D), the particle motion becomes conservative 
significantly simplifying the problem. By demanding V{A) = we find the relationship 
A = ln(l + A^)/Al^ — 1/(1 + A^). It can be established analytically that the growing slope 
of this relation defines the amplitude of the soliton Am- Corresponding nonlinear dispersion 
relation A = \{Am) is exhibited in Fig. 2 (curve ID). One can see from these figures that 
the ID solution is restricted from above: Am < Amcr — 1.4506 . For small amplitudes Eq. 
(jH]) reduces to the standard NSE with a cubic nonlinearity, and ID soliton solution can be 
found analytically (see for cold plasma case when n = n{T)). 

Profiles for large amplitude ID solitons are exhibited in Fig.3. One can see, that as 
A Amcr (A — > A^J,'^^) the profile of the central part of the sohton flattens and widens 
at the top. The existence of fiat-top sohton can be explained by the peculiarities of our 
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+ V{A) 
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focusing-defocusing nonlinearity: the top part of the pulse with A > 1 lies in the defocusing 
region with a tendency for diffraction while the wings of the soliton are in the focusing 
region preventing the total spread of the pulse. It is interesting to remark that for A > A^J.^-* 
the system supports existence of the dark soliton which is an antisymmetric function of 
coordinate with zero intensity at its center. The dark soliton corresponds to the particle 
starting at the right maximum of the curve "c" and going asymptotically towards the left 
maximum of the potential. Background intensity of the field l^ol is bounded from below 
(I^IqI > Aracr) but is uot restricted from above. Thus, a dark soliton with arbitrarily strong 
background intensity is possible in ID. At the critical frequency shift (A = A^^^-*) the dark 
and bright solitons can coexist. The coexistence of these solitons is mainly due to the 
particular type of nonlinearity encountered in the present model. Detailed discussion of the 
properties of dark solitons is beyond the intended scope of this paper. 

In 2D and 3D, the nonzero "friction force" forces the corresponding critical values of 
A to be less than A^^"^-* found for the ID case. In Fig. 2, the relevant curves correspond 
to numerically obtained dispersion relations for 2D and 3D, are displayed. Fundamental 
solitary solutions (without zero nodes) are shown in Figures 4 and 5; In all these examples 
one fundamental feature of the soliton persists, namely that near the critical eigenvalues, 
the profile is endowed with the fiat-top shape. 

The stability of the obtained solutions can be tested by applying the Vakhitov & 
Kolokolov criterion (see 3l| and references therein) according to which the solition is stable 
if dN/dX > , where N = J dr±dC,\A\'^ is the soliton energy ("photon number"). We 
found that in ID the "photon number" is always a growing function of A implying that such 
solitons are stable against small perturbations. 

In Fig. 6 we show the dependence of the "photon number" on the amplitude A^ for 2D 
and 3D solitons. In contrast to the ID case, for the soliton to exist in higher dimensions, its 
energy must exceed a certain critical value N^r- For a 2D case, the "photon number" 
must exceed the threshold energy Ncr = 11.6 for the EM wave to enter the self-guiding 
solitonic regime. Since dN/ dAm > , then dN/ d\ > because Am is a growing function 
of A (see Fig. 2) . Thus the 2D soliton with its power above the critical power is always 
stable. In 3D dN/dAm > provided N > N^r = 236.8 and A^ > 0.6. 



We have demonstrated that the pair plasmas with "assymetry" in initial temperatures 
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of its constituents can support stable large amplitude localized EM wave structures. These 
structures, available in arbitrary (1-3) dimensions, have fiat-top shapes for certain range of 
parameters. This result is particularly interesting for laboratory conditions. Since even a 
"small" difference (very small!) in temperature may be as effective in symmetry breaking as 
the Baryonic correction in early universe, one could readily engineer the laboratory plasmas 
to mimic the cosmic conditions. 

According to fl34l) the total plasma density in the pulse localization area increases dra- 
matically for relativistically intense pulses. Notice, that the general results of the analysis, 
as well as Eq. (jUl) are valid both for ultra-relativistic and non-relativistic temperatures, 
and hence warrants applications to both astrophysical and laboratory plasmas. 



V. MASS ASYMMETRY BETWEEN SPECIES 



In this sub-section we explore another obvious source of asymmetry between species - a 
slight difference in the masses of positive and negative-charged particles. This mechanism 
is different from the one discussed above; the mass asymmetry is initially given and is fixed 
as distinct from the dynamical asymmetry created due to temperature differences. Such a 
plasma can be created by the injection of appropriate ion-beams into a trap. Electron-hole 
plasma in certain semiconductors or e-p collider plasma created by slightly different Lorentz 
factor beams are also possible examples of such system. Appropriate conditions for such 
plasma production could readily appear in Dusty plasmas as well as in astrophysical Jets, 
and Pulsar magnetospheres. 

Much of the framework for investigating this mechanism has already been described. For 
an unmagnetized plasma, one can use the dimensionless equations ffTOj) and ffTTl) in which 
the Lorentz factors of negatively and positively charged particles are, respectively, changed 



1/2 



1/2 



to 7^ = 1 + (p^)^ ' and 7"'' = 6^ + (p^)^ ' , with 6 = /m~ being the ratio of 
positive (m"*") and negative (m~) particle masses. 

For a circularly polarized EM wave with wave frequency cJq = /eq + (1 + 5"^) one can 
readily show that = =i=A (the equivalent of f|T6|) ). Then the standard treatment for 
a transparent plasma (cuq ^ (1 + S^^Y^'^, Vg = ko/ujQ ~ 1) yields the following relations 
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(equivalent of (|25|), (|26|) and (|2 



and 



leading to 



n 



n 

1 + 



1 



l + <? 

l + \A\ 
(1 + 0) 



rr 
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1 + 
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(47) 



(48) 



2iuj— + ^ 



-1^ 



1 



+ VI A+ 



0. 



(49) 



1 + (5((5-0)_ 

For small mass differences 5 = 1 + r] [r] <^ 1) and for the characteristic length L ^ 1, the 
potential may be calculated to be 

explicitly displaying that is proportional to r], i.e., <^ 1 for r] <^ 1 (compare with Eqs. 
f l33|) and fHT]) ). Equations fj49l) with fl50l) . with appropriate normalization of the variables 
and inclusion of transverse field variations, constitute an NLSE (H4l) with the following 
saturating nonlinearity function [r] will be absorbed in the normalization): 

1 



F{\A\ 



1 - 



I 2^ 2 



(51) 



Notice that this form of saturation nonlinearity function coincides with the one obtained in 
24 1 for e-p plasma with a small fraction of heavy ions. The reader may consult Ref. 25I for 



detailed analysis of the system; here we will just summarize the salient features. Equation 
(jHj) (with (I5T1) ) admits a "spherically" symmetric solitary wave solution, i.e. the "light 
bullet", a concentration of mass and energy. And if the "bullet" Am > 0.7, this "bullet" 



is stable. These "light bullets" are found to be exceptionally robust [3^: they can emerge 
from a large variety of initial field distributions and are remarkably stable. The total plasma 
density variation associated with the soliton 6n ~ A"^ is large for A^ ^ 1 ; The solitons with 
ultrarelativistic amplitudes create a large concentration of density. 

Thus, the system of pair plasma with slight initial mass asymmetry between species 
supports the existence of long lived objects - light bullets which carry large amounts of 
mass and energy exactly the same way as the pair plasmas with small fraction of heavier 



10ns 



24M 



26|. 
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The saturating nonlinearity (H5!) caused by an initial temperature asymmetry seems to 
be of a new type (vanishing for intense pulses); it differs quite fundamentally from the one 
found in 2^ (identical to the one originating in a initial mass asymmetry). One would 
think that the temperature asymmetry, manifesting finally as "effective mass" asymmetry, 
would be qualitatively similar to the initial mass asymmetry. But it is not so. Fortunately 
we could trace the cause of the difference. 

The first step in the chain does appear to translate the temperature asymmetry (see flMl) ) 
into a difference in the "effective masses" G^m^". But in reality this "mass asymmetry" is 
dynamical and dependent on the scalar potential 0, rn'^Q/m'^Q ~ (1 — e + O.5z/0) while the 
one originating in an initial mass asymmetry (5 = 1 + r/) is constant in space-time and never 
leads to heating/ cooling. It is of utmost significance to realize that the electrostatic potential 
(f) is important not only for maintaining the nonlinearity, but also to create dynamical 
temperature asymmetry for hot plasma conditions; it is the latter consequence that may 
lead to qualitatively new and interesting phenomena in such a state of matter. 

As a general conclusion we must state that the results of Sec.V are valid only for cold 
plasmas. For initially hot plasmas one must take into account the temperature inhomogene- 
ity effects (discussed earlier in the paper) that change the nature of the nonlinearity and 
hence the localization characteristics of the EM waves. However it is worth mentioning, 
that the main property of density bunching and energy localization is always there in pair 
plasmas with different type initial asymmetries; it is just the character of localization that 
changes with origin of the asymmetry. 



VI. SUMMARY AND CONCLUSIONS 



The main result of this paper is that a fundamentally new type of saturating focusing- 
defocusing nonlinearity is derived for a physical system. This composite nonlinearity, origi- 
nating in a small temperature asymmetry in the constituent fluids of a pair plasma, promises 
the existence of interesting structures that intense electromagnetic waves can acquire in such 
plasmas. We have discussed an illustrative example where different parts of a high amplitude 
pulse are effected differently - the simultaneous expansion of the peak region and scrunching 
of the wings imparts a fiat top shape to the pulse. Most consequences of this new type of 
nonlinearity are yet to be worked out. 
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FIG. 1: The "effective potential" versus the amplitude for different values of the nonlinear 

frequency shift A. The curve "a" corresponds to A > Air^'* ~ 0.2162, the curve "b" has A = X^cr^\ 
and for the curve "c" < A < aIJ.^^ . 
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FIG. 2: Nonlinear Dispersion relations: the effective eigenvalue A as a function of Ajn- The 
boundary line (dotted) corresponds to critical value A = A^J-^^ analytically found only for ID. The 
other three lines represent, respectively the ID, 2D and 3D dispersion relations. 
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FIG. 3: Stationary soliton solution for ID for different critical eigenvalues. Plot "a" corresponds to 
Act- = 0.19315 with Am = 1; plot "b" corresponds to Acr = 0.21583 with Am = 1.4 and plot "c" 
corresponds to Acr = 0.21622 with A^, = 1.47 respectively. The plot "c" represents the flat-top 
soliton solution. 




FIG. 4: Stationary soliton solution for 2D for different critical eigenvalues. Plot "a" corresponds 
to Acr = 0.12789382 with A^ = I; plot "b" corresponds to Acr = 0.17891793 with Am = 1.4: and 
plot "c" corresponds to Acr = 0.20299496 with Am = 1.57 respectively. The plot "c" represents 
the flat-top soliton solution. 
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FIG. 5: Stationary soliton solution for 3D for different critical eigenvalues. Plot "a" corresponds 
to Xcr = 06774722 with Am = 1; plot "b" corresponds to Xcr = 0.12451945 with = 1.4 and 
plot "c" corresponds to Xcr = 0.19222242 with = 1.67 respectively. The plot "c" represents 
the flat-top soliton solution. 
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FIG. 6: The dependence of the "photon number" N on the amplitude A^, for 2D and 3D. Nor- 
malized "photon numbers" correspond to iV^^^) (solid line) and 10~^Ar(^^) (dashed line). For 
2D (3D) the threshold energy for the existence of soliton is Ncr = 11.6 (236.8). 



